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Abstract 

We reply to a recent comment by H. W. Diehl and M. Shpot ( cond\ 



mat/0106501 ) criticizing our paper J. Phys. A: Math. Gen. 34 (2001) 
L327-332. We show that the approximation we use for evaluating higher-loop 
integrals is consistent with homogeneity. A new renormalization group ap- 
proach is presented in order to compare the two methods with high-precision 
numerical data concerning the uniaxial case. We stress that the isotropic 
behaviour cannot be obtained from the anisotropic one. 
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In a recent paper Diehl and Shpot (DS) criticized a method proposed earlier in to 
calculate the critical exponents u L2 , f] L2 , and Jl at order 0(e|) for systems presenting an 
m— fold Lifshitz point. Working entirely in momentum space we perform the calculations by 
using normalization conditions along with dimensional regularization. The symmetry point, 
used to define the renormalized theory, was chosen by setting the external momenta scale 
along the quartic (competing) directions equal to zero, while keeping the external momenta 
scale along the quadratic directions. A detailed account of this approach was given in the 
second paper of reference || for the uniaxial (m = 1) case. 

In momentum space the most general solution to Feynman integrals involving quadratic 
and quartic external momenta scale subspaces perpendicular to each other is a difficult task, 
even at the one-loop level. Indeed, the one-loop integral I 2 contributing to the coupling con- 
stant can be performed exactly only if the external quartic momenta is set to zero. Keeping 
both external momenta scales, one can solve the integral over the quadratic momenta as a 
function of the external quadratic momenta, by choosing Schwinger parameters, for instance. 
The integral over the quartic momenta cannot be obtained as a function of the quartic ex- 
ternal momenta in a closed form. Setting the quartic external momenta to zero simplifies 
the problem, since this integral will contribute with a simple factor to the remaining para- 
metric integrals, which can be solved in a straightforward manner. Absorbing a convenient 
geometric angular factor, the result can be cast in a form which resembles the ordinary (p A 
theory, with replacing e and a slightly different coefficient for the regular term in e^. The 
result is a homogeneous function of the quadratic external momenta scale. 

The parametric integrals play a interesting role in our approximation. To see this, con- 
sider the simplest two-loop integral contributing to the two-point function, namely h(p, k') 
given by 

( k , )= [ d d ' m q 1 d d - m q 2 d m k 1 d m k 2 
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Setting k! = 0, the integral can be evaluated as outlined in |2|]. Before making our 
approximation, one can choose to integrate first either over the loop momenta (qi,ki) or 



over (q 2 , k 2 ). The loop integrals to be integrated first are referred to as the internal bubbles. 
By solving the integral over q 2 first, we obtain 

f°° r d ai da 2 ( ai + a 2 )^ 1 exp( ^-( qi + p) 2 ) [ d m k 2 e~ a ^ ? e~ a ^ k ^ k ^ . (2) 
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Now we can consider the approximation. In order to integrate over k 2 , we have to 
expand the argument of the last exponential. This will produce a complicated function of 
a±,a 2 ,ki and k 2 . Unfortunately, this function has no elementary primitive. Considering 
the remaining terms as a damping factor to the integrand, the maximum of the integrand 
will be either at ki — or at k x = —2k 2 . (The most general choice k\ = —ak 2 yields a 
hypergeometric function.) The choice k\ = —2k 2 implies that k\ varies in the internal bubble, 
but not arbitrarily. Its variation, however, is dominated by k 2 through this constraint, which 
eliminates the dependence on k\ in the internal bubble. At these values of ki, the integration 
over k 2 produces a simple factor to the parametric integral proportional to (cui + a 2 )~~ . 
This allows one to perform the remaining parametric integrals in a simple way. Thus, 
the constraint is designed to preserve the form of the parametric integrals. After realizing 
these integrals, they produce the factor ((qi + p) 2 )~^. Note that the diagrams I 3 and 
I 5 contributing to the two-point function receive the factor after integrating over the 
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quadratic momenta in the external bubble. This factor will not be present in the isotropic 
case, since there is no integration over quadratic momenta to be done in this case. The 
resulting solution to h(p, 0) is a homogeneous function of the external momenta p, not a 
generalized homogeneous function, given by: 

8 — m €l 

The implementation of this constraint on higher-loop integrals proceeds analogously. 
The constraint turns all these integrals into homogeneous functions of the external quadratic 
momenta scale. One can then choose the symmetry point as p 2 = k\ , for example, in order to 
define the renormalized vertices via normalization conditions. The normalization constants 
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Z^(ki), Z$2.{k\) and the beta function are defined in [Q], which give origin to the exponents 
VL2 and r\ui (under the momentum flow in the scale Ki), along with all the scaling relations 
relative to exponents perpendicular to the competing axes ||. This follows in complete 
analogy to the usual (ft 4 theory describing the Ising model. As the constraint is based on a 
physical principle (homogeneity), we do not agree that the approximation is unacceptable. 

Actually, we agree that the speculation made in [§] ( "... suggests that calculations 
performed in momentum space and coordinate space are inequivalent, as far as the Lifshitz 
point is concerned") was unfortunate. (DS) work was the first one to extend the treatment 
done for the m = 2, 6 cases to general m by making use of the scaling form of the free 
propagator in coordinate space representation. However, this does not give them support for 
their speculation that " there is no way that AL's and our calculation can be both correct". 
In the following we shall outline a new renormalization approach in momentum space of this 
problem using two independent fixed points, which is different from (DS) treatment using 
only one fixed point and will help the subsequent discussion. 

So far, we have obtained half the solution to the problem, as we still have to devise a 
method to calculate the critical exponents along the competing axes. We follow a suggestion 
made by Wilson in the early seventies || in order to obtain these exponents independently 
from the ones perpendicular to the competition axes. We can consider another independent 
set of normalization conditions defined at zero quadratic external momenta and nonvanish- 
ing quartic external momenta scale Ki 0. At the Lifshitz point, the free propagator have 
only quartic momenta along the competition axes. Thus, it is possible to perform a di- 
mensional redefinition of the m-dimensional subspace by considering the associated quartic 
momenta to have half the dimension of a conventional momentum scale. As a consequence, 
the term in the bare Lagrangian which is proportional to the quartic momenta does not 
need to be multiplied by another dimensionful normalization constant ((To) in order to be 
meaningful on dimensional grounds. Under a flow in K2 at the corresponding fixed point, 
the normalization constants Z ^(k^), Z 02 (k 2 ) lead to the critical exponents f]L4,^L4 and new 
scaling laws along the competing axes, which are independent of the ones obtained in the 
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subspace perpendicular to the competition directions. In this case, in order to evaluate loop 
integrals we use approximations which preserve homogeneity of the Feynman integrals in the 
external quartic momenta scale k 2 , such that scaling theory is fulfilled. Specifically, consider 
the one-loop integral J 2 (0, P): 



The simplest approximation for this integral which preserves homogeneity in the quartic 
external momenta scale is {{k + P) 2 ) 2 = ((A;) 2 ) 2 + ((P) 2 ) 2 . Of course, more involved approx- 
imations can be developed which preserve homogeneity, but we concentrate on this one for 
the sake of simplicity. The result for this integral is I 2 (0, P) = (((P) 2 ) 2 )~^^-(l + [U] m eL), 



independent infrared divergence of this integral on the external momenta scale along the 
competition axes. The beta function for this case, f3(u) = — 2e L (^f-) K2 , is different (and 
independent) from the one associated to critical exponents perpendicular to the competition 
axes, even though both have the same fixed point at one-loop level. It can be easily checked 
that at the one-loop v Li = ^r. Thus, homogeneity is the guiding principle for obtaining the 
solution to arbitrary loop integrals as a function of k 2 . The resulting scaling relations for 
exponents associated to correlations perpendicular to the competing axes are independent 
from the ones along the competition axes. 

The renormalization group just described can be adapted to treat the isotropic behavior 
m = d close to 8. However, the isotropic case is intrinsically different from this renormaliza- 
tion group perspective. There is only one momenta scale k 2 and just one set of normalization 
conditions. The beta function (3{u) = —£l(^§^)k 2 is na ff the value of the one associated to 
the k 2 characterizing the competing directions in the anisotropic case. They are different, 
since the coupling constant in both cases have different canonical dimension. Technically, 
the isotropic loop integrals do not receive contributions from the parametric integration 
over the quadratic momenta subspace, for they are absent in this case. That is why the 




(4) 



where the geometric angular factor |r(^)r(2 — ^ ; )S' rf _ m S' m has been absorbed in a redefi- 
nition of the coupling constant, and [z 4 ] m = |[1 + ^(1) — tp(2 — ^)]. This result reflects the 
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results for the anisotropic behavior described in 0] cannot be extended to the isotropic one. 
The isotropic behavior has its own scaling relations, which are independent from the ones 
concerning the correlations along the competition axes for the anisotropic case [|3] . 

We can now analyse the previous RG formalisms possessing only one independent mo- 
menta scale for the anisotropic case. The first modern treatment in terms of 1PI vertex 
parts was given by Mergulhao and Carneiro ||. There they set up the formalism in terms 
of normalization conditions in momentum space. They chose the symmetry point at non- 
vanishing quartic external momenta and zero quadratic external momenta as well as two 
conditions on the derivative of the two-point function at two independent external momenta 
scales. This reproduces the earlier scaling relations derived by Hornreich et al |J. They 
went to coordinate space in order to calculate the exponents for the cases m = 2,6. The 
novel feature of this approach is the introduction of an additional normalization constant o"o, 
neeeded to obtain the exponents ^l4,??L4, etc. In 0, (DS) followed this treatment entirely 
in coordinate space in order to extend the formalism to the general m-fold behavior. They 
introduced another normalization constant po m order to treat the crossover and identified 
the critical exponents using the renormalization group in coordinate space. The semianalyt- 
ical coefficients in the e^-expansion are integrals (generalized homogeneous functions) to be 
performed numerically in coordinate space. These numerical integrals only make sense if one 
splits the integration limits on the variable v = a ^xj_ using the scaling and related func- 
tions in the coordinate space representation in the integrand up to the maximum value of 
|v| at I v | = 9.3, and replacing the asymptotic values of these functions for greater values of 
v |7| . This numerical approximation is needed in order to obtain reasonable numerical values 
for the exponents r]L 2 and 77^4 at 0(e|). Otherwise, the generalized homogeneous functions 
in the integrand are not suitable to describe properly the coefficients of the e^-expansion. 
They calculated I2 and I3 along these lines in reference J7|. 

After that, (DS) went to momentum space in order to calculate the two-loop integral 
J 4 (Q,K) using dimensional regularization along with minimal subtraction ||. They used a 
mixed treatment, calculating some integrals in momentum space, going to coordinate space 
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whenever it was convenient (and vice-versa) and making use of a scaling function $ (defined 
in equation (8) of 0). The integral I 2 is a subdiagram of J 4 , depending on two external 
momenta scales as well. Nevertheless, they fixed the quartic external momenta scale to be 
zero and concluded that I4 does not depend on it for general values of the quartic external 
momenta scale. Indeed, according to equation (B.14) of |5J (e^ = 4 + — — d) 



h{Qe±, K) = 9-1 [- + J u (m) + 0(e L )\ . (5) 



This happens to be incomplete. The problem can be traced back to the calculation of the 
one-loop integral: 

d d - m qd m k 



I 2 (P±,K(,) = 7 - ^ r . (6) 

" J {(k + K')l+(q + P)l)(kl + ql) 



They only computed this integral for vanishing external momenta along the quartic direction 
(Ki =0). In this case, one has 

I 2 (P ± ,K' ll =0) = (P 2 )-^I 2 (e ± ). (7) 

However, by setting P± = 0, keeping K'» different from zero, we have from the discussion 
following Eq. (4) 

I 2 (P ± = 0,K'e l] ) = (K' i )-^h(e ll ), (8) 

We stress that /2( e ||) and /2( e ±) are different functions in general. One can choose them to 
have the same leading singularities (the multiplicative factor is absorbed in a redefinition of 
the coupling constant anyway), the difference appearing in the regular terms in e^. In fact, 
the complete integral I 2 depends on these two momenta scales. 

We recall that in a proper minimal subtraction procedure, all the external momenta 
should be kept arbitrary 0. (DS) did not take into account this fact to proceed with the 
minimal subtraction, rather keeping only the quadratic external momenta Qe± in I4 and 
P± in I 2 , and setting the quartic external momenta to zero in these integrals. They should 
show how the necessary cancellations of poles take place along the quartic subspace as well 



in order to have a satisfactory minimal subtraction scheme. In fact, the cancellations along 
the quadratic directions actually work when one fixes the quartic external momenta to zero. 
Even though this procedure is not complete, one can accept it as a new type of minimal 
subtraction to this problem. 

For the anisotropic case, the (DS) method is based on one fixed point and the (almost 
exact) numerical integration for two-loop integrals which appear as the coefficients of the 
e^-expansion using coordinate space representations whenever it is convenient. On the other 
hand, the method developed in reference || (and discussed here) in momentum space utilizes 
two fixed points. This new method states that there are four independent critical exponents 
(instead of three) with two independent set of scaling laws relating exponents along the 
quadratic and quartic directions in each subspace separately. We use an approximation, 
namely the constraint relating loop momenta in internal and external subdiagrams, which 
yields analytical results for higher-loop integrals. 

A comparison of the two methods with numerical results for the exponents associated to 
perpendicular correlations to the competing axes (labeled with the subscript L2 after 
is in order. For the m = 1, d = 3, JV = 1 case, (DS) found using MATHEMATICA: 
v li — 0.71,7^2 = 1-40. This is consistent with the newest Monte Carlo simulations for 
7L2 = 1.36 ±0.03 [f[C)|l , and compatible with an earlier Monte Carlo study (1.40 ±0.06) [jTTjl . 



On the other hand, our approximation yielded v L2 = 0.73 and 7^ 2 = 1-45 @. When 
using the new hyperscaling relation obtained in || for the specific heat exponent, namely 
2 — a L2 — (d— y) u L2 and replacing the value v L2 = 0.73, we obtain a i2 = 0.175, whilst the 



most recent Monte Carlo calculation is atm = 0.18 ±0.02 |T(|. On the other hand, using the 
value obtained by (DS) vl2 = 0.71 in the new hyperscaling relation, we find ai2 = 0.225. We 
can proceed and analyse the new scaling law obtained in || for the magnetization exponent 
Pl2 = 2 PL2 ((d ~~ Y ) ~~ 2 ± f]L2)- Our calculation yields (3l2 = 0.198, whereas the simulation 
result is (3l2 = 0.238 ± 0.005. Using (DS) results for ui2 and r]L 2 inside this new scaling 



relation one finds (3l 2 = 0.192. The high precision numerical values |TD[ are in very good 
agreement with the two-loop results using our approximation, which we believe cannot be 
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said to be "unacceptable" at this point. The very similar values obtained for the exponents 
using either (DS) or our two-loop calculations confirms that momentum and coordinate space 
calculations should give the same results, since either our approximation or (DS) numerical 
approximation is responsible for a rather small deviation in the two results when compared 
to the above numerical output. 

Finally, we emphasize the failure of (DS) method to treat the isotropic behavior. As was 
pointed out in and explained in this work, UL2 and r]L2 are not valid for the isotropic 
(m = 8) case. At the Lifshitz point, only the momenta scale along the competing axes is 
meaningful for the isotropic case. Hence one has to start from scratch using this momenta 
scale, which is incompatible with our choice of normalization conditions. (DS) calculated 
I 2 and I4 along the components of quadratic external momenta only. But this momenta 
scale makes no sense for the m = 8 case at the Lifshitz point, since they are not present any 
longer. In that case the fixed point should be determined entirely as a function of the quartic 
external momenta scale as shown in ||. There it was found that the isotropic behaviour 
cannot be obtained from the anisotropic behaviour. This is in contradiction to (DS) and we 
conclude that it is most likely the use of the momentum scale vanishing at m = 8 that led 
(DS) to erroneous results. 

In conclusion, we have shown that our two-loop results do constitute a very good ap- 
proximation for calculating critical indices. Our method proved to be very simple to give 
analytical expressions to the exponents. It is based on a renormalization group analysis 
consisting of two independent fixed points and is a natural alternative to the (DS) semian- 
alytical approach based on only one fixed point. In view of the comparison with numerical 
values, we believe that both methods for the anisotropic cases deserve further investigation 
in order to unravel the fascinating issues concerning the Lifshitz critical behavior. 

The authors would like to thank B. V. Carlson for a critical reading of the manuscript 
and support from FAPESP, grant numbers 00/03277-3(LCA) and 00/06572-6(MML). 
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